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Abstract 

We study the evolution of linear density perturbations in a large spherical void universe which 
accounts for the acceleration of the cosmic volume expansion without introducing dark energy. 
The density contrast of this void is not large within the light cone of an observer at the center of 
the void. Therefore, we describe the void structure as a perturbation with a dimensionless small 
parameter k in a homogeneous and isotropic universe within the region observable for the observer, 
j We introduce additional anisotropic perturbations with a dimensionless small parameter e, whose 

evolution is of interest. Then, we solve perturbation equations up to order Ke by applying second- 
order perturbation theory in the homogeneous and isotropic universe model. By this method, we 
can know the evolution of anisotropic perturbations affected by the void structure. We show that 
the growth rate of the anisotropic density perturbations in the large void universe is significantly 
different from that in the homogeneous and isotropic universe. This result suggests that the 
observation of the distribution of galaxies may give a strong constraint on the large void universe 
model. 
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I. INTRODUCTION 



In standard cosmology, the Copernican principle is applied i.e., it is assumed that we are 
not living at a privileged position in the universe. Combining the Copernican principle with 
the isotropy of the Cosmic Microwave Background (CMB) radiation, leads to the conclusion 
that our universe is well described by the homogeneous and isotropic universe model. In 
the framework of the homogeneous and isotropic universe, the observational data of the 
luminosity distances of Type la supernovae (SNIa) indicates an acceleration of the cosmic 
volume expansion of our universe. The acceleration of the cosmic volume expansion in the 
homogeneous and isotropic universe implies the existence of so-called dark energy that acts 
as a source of a repulsive gravitational force if we assume general relativity at cosmological 
scales. At present, there is no theory that can naturally explain the origin of dark energy, 
and it seems worth investigating alternative scenarios. In order to do this, we have to discard 
general relativity or the homogeneity assumption. 

Inhomogeneous cosmological models without dark energy have been proposed indepen- 
dently by Tomita and Celerier Q. In Celerier's model, the observer is located at the 
symmetry center of a very large spherical void which can explain the SNIa observations. 
Since the observer is located at a special position in the universe, we call this model a "the 
non-Copernican universe model" in this paper. The common assumption of non-Copernican 
universe models is that an observer is located in the vicinity of the symmetry center, which 
explains the fact that the CMB radiation is observed to be isotropic. The most common way 
to describe non-Copernican universe models is to use the Lemaitre-Tolman-Bondi(LTB) so- 
lution for the Einstein equations, which describes the motion of spherically symmetric dust. 
SNIa observations in non-Copernican universe models have been studied by many re- 



searchers [l 
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^ and it has been proven that the distance- redshift relation in the ACDM 

universe can be reproduced using LTB universe models , 0, Esf. Non-Copernican uni 



verse models have been tested by other observations including the CMB acoustic peaks 



24J |. the radial baryon acoustic oscillation scales |24H26|. the kinematic Suniyaev-Zeldovich 



effect 27H31| and others 32N46| . Although these observations have imposed restrictions on 
these models, they have not yet ruled out the models. It is not easy to confirm whether our 
universe follows the Copernican principle. 

In this paper, we focus on the evolution of structures such as clusters of galaxies and super- 
clusters in the non-Copernican universe in the matter dominant era which is well described 
by the LTB solution. It is expected that observations of the large-scale structures and their 
evolution can be used to test the non-Copernican universe model, since the evolution of the 
anisotropic perturbations reflects the tidal force in the background spacetime. However, the 
evolution of perturbations in the LTB solution has not yet been fully studied. This is because 
the isometrics in the LTB spacetime are less than in the homogeneous and isotropic universe. 
Although master equations for perturbations for general spherically symmetric spacetimes 
have been derived a long time ago 47| (see also Ref. 48| for the LTB background), these 
equations for the LTB solution cannot be reduced to ordinary differential equations, in 
general. This is a very different situation from the case of the homogeneous and isotropic 
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universe. 

Recently, Alonso et al. [49[ performed numerical simulations for non-Copernican models 



including only cold dark matter. They studied the perturbed Einstein-deSitter universe 
with two kinds of perturbations: one forms a spherical void, and the other is a non-spherical 
perturbation with a power spectrum with a random phase Gaussian probability distribu- 
tion. They followed the growth of these perturbations using Newtonian A^-body simulations. 
However, in order to confirm the validity of the tnumerical simulations, analytic comple- 
mentary studies are necessary. Some authors 2l|, |35|, l50|] have studied perturbations by 
using a "silent approximation" that neglects the magnetic part of the Weyl tensor. Here, 
we should note that the magnetic part of the Weyl tensor usually plays an important role 



even in Newtonian situations 5l|- Hence, in this paper, we propose another complementary 
analytic approach. 

In many non-Copernican models, the void structure becomes nonlinear at the present 
time. However, Enqvist et al. [sgI pointed out that the void inhomogeneity remains in a 
quasi-linear regime ~ (9(0.1) inside a past light-cone of an observer at the center of the 
void. Actually, they considered a linear perturbation in the Einstein-deSitter universe that 
is consistent with the SNIa data, and showed that the fraction of the spherically symmetric 
linear perturbation does not exceed 30% inside the past light-cone. This result implies that 
non-Copernican LTB cosmological models compatible with the observed distance-redshift 
relation may be studied by perturbation theory for the homogeneous and isotropic universe 
filled with dust at least for the inside of the past light-cone of the central observer. 

In this paper, we investigate the growth of perturbations in the non-Copernican uni- 
verse models by applying the above idea. It is rather difficult to analyze the evolution 
of anisotropic perturbations in the non-Copernican LTB universe model, while it is much 
easier to study the evolution of non-linear perturbations in the homogeneous and isotropic 
universe model by successive approximation. We adopt the latter approach. We introduce 
two-parameter perturbations with small expansion parameters k and e in a homogeneous 
and isotropic dust universe. The limit e — ?■ leads to the exact LTB solution, if we take 
all orders of k into account. By contrast, the limit k — )■ with < e ^ 1 leads to the 
homogeneous and isotropic universe with small anisotropic perturbations. Then, in order 
to see the effect of the void structure on the evolution of the anisotropic perturbations, we 



study the non-linear effects up to the order of Ke, following Ref. 52 



This paper is organized as follows. In § HIl we derive the equations for perturbations 
parametrized by k and e in the homogeneous and isotropic dust universe and obtain gen- 
eral solutions up to order ne. In § llllt by fixing the initial conditions, we calculate the 
angular power spectrum of the density perturbations. In § IIV[ we analyze the growth of the 
perturbations by using the angular growth rate. §|V]is devoted to a summary and discussion. 

In this paper, we use the geometrized units in which the speed of light and Newton's grav- 
itational constant are one, respectively. The Latin indices denote the spatial components, 
whereas the Greek indices represent the spacetime components. 
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II. TWO-PARAMETER PERTURBATIONS IN A HOMOGENEOUS AND 
ISOTROPIC DUST UNIVERSE 



A. Perturbations with two kinds of parameters 

As mentioned, we study the perturbations in the homogeneous and isotropic universe 
which is often called the Friedmann-Lemaitre- Robertson- Walker (FLRW) universe. Since the 
structure formation begins after the universe has begun to be dominated by non-relativistic 
matter, it is sufficient for our purpose to consider the universe model filled with dust. Using 
the spherical polar coordinates for 3-dimensional space, the line-element is given by 

= -dt^ + a^{t) (^Y^-j^ + r^dn^^ =: -dt'^ + a\t)-fijdx'dx^ , (1) 

where a{t) is the scale factor which will be determined by the Einstein equations, K is 
constant, dQ"^ is the 2-dimensional round metric, and, for later convenience, we have defined 
the background conformal 3-metric jij. The constant K has the same sign as that of the 
curvature of the 3-dimensional space specified by t =constant. The stress-energy tensor of 
dust is given by 

f"" = p{t)u^u'', (2) 

where p{t) is the energy density, and is the 4-velocity whose components are given by 
= (1, 0, 0, 0). The Einstein equations for the FLRW universe are 

d\ Sir K d Air _ 

= -irP - ^ and - = — -p, (3) 



a 

where a dot denotes differentiation with respect to t. 

Since our main interest is the evolution of density contrasts and their correlations, we 
consider only scalar perturbations in the FLRW universe. As mentioned in the previous 
section, we introduce two small independent non-negative parameters, k and e. The limit 
e — )■ leads to the exact LTB solution, if we take all the orders of k into account. By 
contrast, the limit k — > with < e <^ 1 leads to the homogeneous and isotropic universe 
with small anisotropic perturbations. 

Then, by choosing the synchronous comoving gauge, the line element of the perturbed 
spacetime can be written in the form 



ds^ = -dt^ + a^jt) K 



N 

N=0 



+ 0{e 



2\ 



(4) 



where /j"^ = /f = 1, r2 = {e,(f)) are the polar and azimuthal angles, and Vi denotes the 
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covariant derivative with respect to 'jij. The perturbed stress-energy tensor is given by 

N=0 

where A'-^^ = 1. If we wish to study the evolution of the perturbed FLRW universe with 
the same accuracy as the hnear perturbation analysis for the LTB solution, we should take 
all orders of k, and the first order with respect to e. However, if k is much smaller than 
unity, it will be possible to evaluate the evolution of the anisotropic perturbations in the 
LTB solution by studying up to the first order with respect to n. In this approximation, the 
effect of the void structure on the evolution of anisotropic perturbations appears at order 



B. First order perturbations 



We expand the perturbation variables in terms of the spherical harmonic functions 
Yimi^)- The Einstein equations of order k correspond to the equations for the perturbations 
of £ = mode: 



(rl 



/(I) 



21 



(1) 



-2A«, 



where 



(6) 
(7) 

(8) 
(9) 



and a dash denotes a partial differentiation with respect to r. The Einstein equations of 
order e lead to the equations for the perturbations oi i > modes, and we obtain 



A'^^{t,r)-KB',:^{t,r 
2(1 - Kr^ 



+ 



r 

3K 



Kr 

+ 1) 



-'im 



0, 



5«(t,r) + 2i75(;^(t,r)-4vrp5«(t,r) 

where we have used the eigenvalue equation 

1 92 



92 cos6' d 

de^ ^ ^^de^ sin2^902 



Y.Jft) 



(10) 

(11) 
(12) 

(13) 



We note that the equations of order e and k decouple with each other. For later convenience, 
we also show the equation for B{t, r, f2) before deriving Eq. (fTTl) by the spherical harmonics 
expansion: 
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(14) 



General solutions for Eqs. ([6]) and (fT2|) are given by 

A(^)(t,r) = D+{t)Al\r) + D-{t)A'^l\r), (15) 

dt,r) = D^t)6fl^{r) + D-it)6fl{r), (16) 

where 

D^{t) = H j and ^-(t) = H, (17) 

and aI^'* and stand for initial values. and D~ represent the growing and decaying 
modes, respectively. 



C. The perturbations of order k€ 

As already mentioned, we are interested in the effect of the void structure on the evolution 
of anisotropic linear perturbation, and this effect first appears at order ne. Hence we shall 
focus on the perturbations of this order. 

The perturbations of order ne correspond to the second order perturbations in the FLRW 
universe model. From the second-order Einstein equations together with the background 
and the linearized Einstein equations, we obtain the evolution equations for the expansion 
coefficients of 5^'^\t,r, ft) with respect to Yim{^) as follows: 

where the source term Sirn{t,r) is given by 

S,™(*. r)= -/(,") (a- - ^) B« + (/« 
+ (2A'« - if) Sg] + 8.pA")i«. 

By solving Eq. ( fTSl) . we obtain 

^ ls,^M (^^M^zfWi^) * (20) 

where W{s) is the Wronskian given by W{s) = D~^{s)D~{s) — D~^{s)D~{s), and homoge- 
neous solutions have been absorbed in 5^^. Then, we obtain the anisotropic linear density 
contrast Sem in the LTB solution as 

6Ut,r) = e5lll{t,r) + Ke 5,„(s,r)(^ ^ ^^^^ ^' —jds 

+ 0{Kh), (21) 

where is the initial time. 



S limit, r), 



(18) 



Im, 



(19) 
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Hereafter, we neglect the decaying modes of order e. By using Eqs. (fTT]) and f|T6|) . is 
written in the form 



B^ll{t.r) = D^{t)BZ{r 



(22) 



where B^l)^{r) is the initial value. Then, the anisotropic density contrast ( 12T]) is rewritten 

as 



>«+/ 



where Ti, T2 and T3 are defined by 

'D-{t)D+{s)-D+{t)D-{s) 



ds 



W{s) 



X 



ds 



W{s) 



2r2 



Ts{t,r) 



ds 



D-{t)D+{s) - D+{t)D'{s) 



xD+{s) {fl\s,r)-l\''{s,r) 



W{s) 



(23) 



D+{s) (2A^^\s,r) -ill\s,r)j +D+{s) x 87rp{s)A^^\s,r^ , (24) 
D-{t)D+{s)-D+{t)D-{sy 



(25) 



(26) 



III. ANGULAR POWER SPECTRUM AND ANGULAR GROWTH RATE 



In the previous section, we derived the growing solutions for density contrasts. Once 
we have specified the isotropic linear perturbations /|*-|^'*, Z^'*, A''^-' and the initial anisotropic 

inhomogeneities S^^^, B^-^, we obtain the density contrasts by Eq. fl23|) . In this section, 
we derive the explicit form of the angular power spectrum of the density perturbation S^rn 
for a given set of initial conditions in terms of the standard power spectrum. Then, we 
define the angular growth rate by using the angular power spectrum. Hereafter, in order to 
determine the perturbations of order k, we refer to the non-Copernican LTB universe model 
with uniform big-bang time (see Appendix A). 



A. Initial power spectrum of density contrast 



By virtue of the uniform big-bang time (see Appendix A), the present non-Copernican 
LTB universe approaches the homogeneous and isotropic universe as time goes back. Hence, 
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it is reasonable to assume that the initial conditions for the anisotropic perturbations are 
the same as in the case of the FLRW universe. Then, the initial power spectrum of the 
density contrast can be expressed as follows: 

k)5(^)(ti, k')) = {2nf6Uk - k')P{U, k), (27) 

where 6d is the Dirac's delta function, ti represents some sufficiently early time already 
introduced in Eq. (12T]) . and the Fourier transform of the density contrast is defined by 

(t, k) = y d^x6^^^ (t, x)e-'^". (28) 

If we choose the initial time after recombination, the matter power spectrum including 
baryons and cold dark matter can be written as 

P{U,k) = [D^{U)]'P{k), 

P{k) = Aok^'T^k), (29) 

where Aq is a positive constant which represents the amplitude for perturbations on large 
scales, n is constant, and T{k) is the matter transfer function. In this paper, we assume 
the Harrison-Zel'dovich spectrum n = 1. As for the transfer function, we adopt the fitting 



formula developed by Eisenstein & Hu |53| (see Appendix B) 



B. Angular power spectrum and angular growth rate 

In order to observationally study the evolution of perturbations in the non-Copernican 
universe, we need to specify the observable quantities by using the density contrast (l23l) . 
The simplest quantity that we can currently calculate is the angular power spectrum. We 
define the angular power spectrum of the density contrast by 



Ci{t,r) 



2i + 



(30) 



m=-e 



where * denotes the complex conjugate. 

Hereafter, we focus on the non-spherical perturbations whose wavelengths are much 
smaller than the spatial curvature radius {k ^ a/I-R'I). Then, from Eq. ( fT^ . we have 

eB^^\t,k) -25^'\t,k). (31) 

The initial values 5^),\(r) and Bf^{r) which appear in Eqs. (l22l) and ( 1231) can be written using 
the Fourier transform of the initial density contrast 5^^\ti,k) as 



\ ^ 1 


(47r/) j 











d^k 



5W(ti,k)j,(fcr)r4(a 



(27r)3 

d^k / 2 

^5«(t,k)j,(A:r)F;„(a) (^^ 



(32) 
(33) 
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where fik denotes the polar and azimuthal angles in the Fourier space, and we have used the 
relation between the expansion coefficient with respect to Y£m{^k) and the Fourier transform 



(34) 



and Eq. ( 13T1) . Using Eqs. (l27j) . (l29l) . (l32l) and (l33l) . the angular power spectrum of the 
density contrast ( l23l) can be rewritten in the following form: 



C£(t,r) = e^D+'^{t)Ki{e,r) + Ke^\2D+{t)Ti{t,r)Ki{i,r) 

+2D+{t)f2{t, r, i)K2{i, r) + 2D+{t)fs{t, r)K,{£, r)] + 0{^'e^] 



(35) 



where 



i^i(£,r) 
i^2(^,r) 







2\ A-oo 







hi: 



dkP{k){krf]j{kr 
dkP{k){krfj1{kr 
dkP{k){krYji{kr)j(^{kr 



(36) 



and T2 and T3 are defined as the short wavelength approximation {k ^ -\/|^) of T2 and T3 
by 



T2(t,r,£) 



(is 



D-(t)D+(g) -Z}+(t)D-(s) 

iv(s) 



xD^{s){i^l\s,r)-i\l\s,r] 



+ 1) 



2r2 



(37) 



T3(t,r) 



ds 



D-{t)D+{s)-D+{t)D-{s) 



W{s) 



xD-^(s){i^l\s,r)-i\}\s,r] 



(38) 



Once the initial density power spectrum P{ti, k) is specified, we can calculate the angular 
power spectrum Ci{t,r) by using Eq. ( 135|) . 

To investigate the growth rates of the perturbations, we define the angular growing factor 

by 



A(t,r) 



1/2 



(39) 



It is easy to see that the angular growing factor Di(t, r) is equal to D^{t)/D^{t-^ up to order 
e. 
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Hi"- 









(41) 



Basically, an observer can see the cosmological structures on his/her past light cone by ob- 
servations through electromagnetic radiation^ and hence, in the case of the non-Copernican 
universe model, it is useful to consider quantities on the light cone of an observer who stays 
at the symmetry center of the void at present. Hereafter, for simplicity, we call the observer 
who stays at the symmetry center of the void at present "the central observer" , and the past 
light cone of the central observer is denoted by Sic. The past light cone Sic is generated 
by the past-directed outgoing radial null geodesies = {dt/d\,dr/d\, 0,0), where A is the 
affine parameter. The cosmological redshift z is defined by 

dt/dX _ 

{dt/dX)o ' ^ ^ 

where dt/dX and {dt/dX)^ are the value at the time of the emission of a photon and that 
at the time of the detection of the photon by the central observer, respectively. By using 
the cosmological redshift z instead of the affine parameter A, the geodesic equations for the 
generator of the past light cone Sic up to order n are given by 

dr Vl — Kr'^ 
dz (1 -|- z)aH 
dt 1 ; \ 

di = -(miff V - m''0 • (^2' 

We denote the solution of the above equations hj t = tydz) and r = ric(z). 

Then, using the angular growing factor, we define the angular growth rate on the past 
light cone Sic as a function of redshift z as follows: 

f f ^ c?[ln A(tic(z),ric(2:))] , . 

^'^'^ dln(l + .) • ^^^^ 

Here, we note that the angular growing factor and the angular growth rate do not depend 
on the amplitude Aq in Eq. ( 129|) . We also note that the angular growth rate up to order 
e agrees with the growth rate usually used in the linear perturbation theory of the FLRW 
universe, d{\nD~^)/d{lna). 



IV. EVOLUTION OF DENSITY PERTURBATIONS IN THE CLARKSON-REGIS 
MODEL 

A. Linearized Clarkson-Regis model 

In order to determine the perturbations of order k, we use the non-Copernican LTB uni- 



verse model given by Clarkson and Regis 17| (see Appendix A), which we call the Clarkson 



^ The past light cone of an observer at the event p is defined by the boundary of the causal past of p, which 
is usually denoted by J~ (p) in general relativity. Strictly speaking, the observer can see the inside of the 
light cone through a congruence of the light rays which have experienced caustics caused by gravitational 
lens effects or scattering due to electromagnetic interactions in the real universe. 
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and Regis (CR) model. We shall study the evolution of linear anisotropic perturbations 
in the CR model by using the second order perturbation theory of the FLRW universe 
filled with dust. In order to approximate the CR model by the linearly perturbed FLRW 
universe filled with dust, we must first specify the background FLRW universe. Here, we 
determine the background FLRW universe so that the cosmological density parameter of the 
background is equal to 0.242, which is equal to the value of the density parameter function 
^m{i^) at the symmetry center of the CR model (see Eq. (]A13P ). 
We define the "density contrast" of the CR model as 



p(CR)(t, 0) 



(44) 



where is the energy density of the CR model. The density contrasts A'^^^^'s on three 
constant time hypersurfaces are depicted in Fig. [1] as functions of r. Here, we have used 



, (CR) 



. 5 




r[Gpc] 



FIG. 1: Density contrasts on the spacelike hypersurfaces for t = tic(lOO), t = tic(l) and t = tQ, as 
functions of r. 

the cosmological redshift z to specify each constant time hypersurface given hj t = tic{z). 
We can see that the void structure grows with time. Since the big-bang time is uniform, 
there is only the growing mode in the CR model. The void size is about 12Gpc, and the 
vicinity of the center is locally the dust filled FLRW model with the cosmological density 
parameter = 0.242, whereas the asymptotic region is almost the same as the dust filled 
FLRW model with = 0.7. The Hubble parameter at the center is Hq = 74kms~^Mpc~^. 

By using the density contrast A^^^\ we give the initial conditions for the isotropic linear 
density contrast A^^ in Eq. ( !T5|) as follows. As mentioned, since the CR model has only 
the growing mode, we should set A^^{r) = 0. By contrast, A^^(r) is determined by the 
assumption that the density contrast A^^^ exactly agrees with that of the CR model at the 
initial time A^^^\ti,r), i.e.. 
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where the initial time is determined by ti = tic(lOOO). Once the initial condition for the 
density contrast is given, we obtain all perturbations of order k, by solving the perturbation 
equations up to the corresponding order, and as a result, the linearized CR model is obtained. 

To evaluate the accuracy of the linear approximation, we plot A^^^^ and A*^^) on the past 
light cone Sic in Fig. |2l The relative error between the exact and linearized CR models is 

. 5 I 1 



0.4- 




FIG. 2: Exact density contrast A'-'^'^^ and that of the linearized CR model A^^^ on the past light 
cone Sic, plotted as functions of the redshift z. 

less than 30% on the light cone Sic. Inside the past light cone Sic, the error is smaller than 
that on the past light cone Sic, since the CR model has only the growing mode. There is 
no qualitative difference between the exact and the linearized CR models, and thus we may 
see the qualitative behavior of linear perturbations in the CR model by the perturbative 
analysis of the FLRW universe based on the linearized CR model. 



B. Evolution of density contrasts in the CR model 



Let us consider the evolution of the anisotropic density contrasts in the CR model. By 
using the angular power spectrum Ce{t, r) given by Eq. (135|) and the transfer function T{k) 
given in appendix we depict the angular growing factors De{t,rys defined by Eq. fl39|) at 
each comoving distance as functions of t in Fig. [31 

Here, we introduce a useful quantity defined by 

~k:=- . (46) 
r 

Note that k is equal to the comoving wave number of the mode £ at a distance r in the flat 



sky approximation (see e.g., Ref. [5J]). In Fig. [3l the present time is Hoto = 0.83, and £ is 
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FIG. 3: Angular growing factors D^'s in the CR model at r = 40Gpc (dotted line), r = 4Gpc 
(dashed line) and r = (dot-dashed line) depicted as functions of t. The present time is //q^o = 
0.83. We choose I so that k = O.SMpc"^ 
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FIG. 4: Angular growing factors D^'s in the dust filled FLRW universe models with f^M = 0.242 
and = 0.7, together with that for the CR model. 

chosen so that k = 0.5Mpc~^. This choice of £ shows us the evolution of perturbations with 
the size of a cluster of galaxies, i.e., 27r/A; ~ lOMpc. 

We can see from Fig. Elthat the larger the comoving distance of a perturbation from the 
symmetry center, the faster the growth of the perturbation. This result may be explained by 
the fact that the energy density of the CR model is a monotonically increasing function of r, 
since the growth rates of perturbations in the FLRW universe is a monotonically increasing 
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function of ^m- We also depict De of the FLRW universe models with = 0.242 and 0.7, 
respectively, together with that of the CR model in Fig. HI We can see from this figure that 
Di of the FLRW universe with Q,u = 0.242 agrees with Di{t,r = 0) of the CR model. We 
note that in the FLRW universe with = 0.7 does not agree with that far from the void 
(r = 40Gpc) in the CR model. This result might not be real, but rather could be an error 
caused by using the linearized CR model, since the error due to the linear approximation 
becomes larger for larger r. 




Hot 



FIG. 5: Angular growing factors Di{t, r = 4Gpc) in the CR model for k = O.OlMpc"^ (solid line), 
O.lMpc^^ (dashed line) and IMpc^^ (dot-dashed line), depicted as functions of t. The right panel 
shows a close-up of the 0.7 < H^t < 0.71 

We have also investigated the dependence of Di on k. The angular growing factors 
Di[t,r = 4Gpc)'s with various values of i, or equivalently, k are depicted as functions of t 
in Fig. O We find that the dependence of D£{t, r = 4Gpc) on k is very small in the case of 
the CR model. 

Next, the angular growth rate fi defined by is plotted as a function of z for the 
CR model, together with those of the dust filled FLRW with = 0.242 and = 0.7 
and the fiat ACDM with Qm = 0.28 in Figs. El and [3 Figure [7] shows a close-up of the 
region < 2; < 1 in Fig. [HI Here, i is also chosen so that k = 0.5Mpc~^. From Fig. El we 
can see that the value of of these models approach the value of the Einstein de-Sitter 
universe (/^ = 1) in the high redshift domain. We can see from Fig. [7] that at the central 
observer 2 = of the CR model agrees with the value of the dust filled FLRW universe 
with Qyi = 0.242. We can also see that in the CR model is significantly different from 
those in homogeneous and isotropic universes for redshift < 2; < 1. This result implies 
that, if we can somehow observe fi, the observational data may give a strong constraint on 
non-Copernican universe models. 
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FIG. 6: Angular growth rate fi in the CR model (solid line), together with those of the dust filled 
FLRW universe with ilu = 0.242 (dot-dashed hue) and J^m = 0.7 (dashed line) and the flat ACDM 
universe with Qm = 0.28 (dotted line) as a function of z. £ is chosen so that k = 0.5Mpc~^. 
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FIG. 7: Close-up of the region < z < 1 in Fig. [6l 



V. SUMMARY AND DISCUSSION 

We have studied the evolution of anisotropic density perturbations in a large spherical 
void model, which is a non-Copernican universe model based on the LTB solution. As is 
well known, the analysis of anisotropic linear perturbations of the LTB solution is much 
harder than that of the homogeneous and isotropic universe. Therefore, we have not tried 
to directly solve linear perturbations in the LTB model. Instead, we have studied linear 
perturbations in a linearized LTB model. The linearized LTB model is obtained by regarding 
the inhomogeneities in the LTB model as isotropic linear perturbations in the dust filled 
FLRW universe. Hence, our method is relevant for the LTB model only if the isotropic 
inhomogeneities is so small that they can be treated as linear perturbations on a dust filled 
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FLRW universe. 

Assuming the uniform big-bang time, we have given the initial conditions for anisotropic 
perturbations, which are the same as those of the FLRW universe model. Then, we calculated 
the angular power spectrum of the anisotropic density contrast taking into account the tidal 
effect of isotropic inhomogeneity, and further we defined the angular growing factor and an 
angular growth rate using the angular power spectrum. 

In § IIVI we calculated the angular growing factor using the linearized Clarkson-Regis 
(CR) model, which has the uniform big-bang time. From the behavior of the angular 
growing factor, we found that the speed of growth of a perturbation is a monotonically 
increasing function of the comoving distance from the center of the void. Because of this 
property, the angular growth rate in the CR model differs from that in the dust filled FLRW 
universe even for low redshift {z < 1). So, if we can observe the angular power spectrum of 
the density contrasts and can observe their angular growth rate, the observational data will 
strongly restrict non-Copernican universe models. 

Guzzo et al. and Blake et al. have presented the observational results on the growth rate 



of the density perturbation as a function of the redshift z on our past light cone 55|, |56j . Here, 



we should note that their observations are a distortion of the correlation function of galaxies 
in redshift space and they derived the growth rate from it assuming the homogeneous and 
isotropic background universe. Hence, our present result regarding the angular growth rate 
does not correspond to their "observational results on the growth rate". To calculate the 
redshift space distortion of the galaxy distribution in non-Copernican universe models, we 
have to consider not only the growth in the transverse direction but also the radial direction. 
This is a future work, and we shall present the results elsewhere. 

We expect that the formula developed here could be applied to many inhomogeneous 
cosmological models, since we have only assumed that the isotropic inhomogeneities are 
small inside the past light cone Eic. We also expect that the formula could be potentially 
used for calculating the integrated Sachs- Wolfe effect, baryon acoustic oscillations, the shape 
of the 3D power spectrum and so on. 
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Appendix A: Clarkson-Regis model 



The Clarkson-Regis (CR) model is a non-Copernican universe model based on the 
Lemaitre-Tolman-Bondi (LTB) spacetimes, which is an exact solution of the Einstein equa- 
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tions and describes the motion of spherically symmetric dust. The CR model can explain 
the SNIa data and peak positions of the fluctuations in the CMB radiation. The line element 
and the stress-energy tensor of the LTB spacetime are given by 



ds 



-de 



{drR{t,r)f 



1 — k{r)r'^ 



■dr^ + R'^{t,r)dn^ 



(Al) 
(A2) 



where k, and p are an arbitrary function of the radial coordinate r, the 4- velocity and the 
energy density of the dust, respectively. The 4-velocity of the dust fluid is given by = 6q. 
The Einstein equations lead to 



R 

p{t,r) 



2M(r) k{r)r'^ 



i?3 

drM{r) 



A-KR^{t,r)drR{t,ry 

where M(r) is an arbitrary function of r. By integrating Eq. flA3l) . we obtain 

R{t,r) = {<6M{r)f'\t - t^{r)f'^S{x), 
where t-Q{r) is an arbitrary function of r, x is defined by 



and, by defining r] as 



X := k{r)r'^ 



-(sinh ^J—r| 



6M(r) 



2/3 



2/3 



X 



Q2/3 



the function S{x) is given by 



62/3 



cosh yf—rj — 1 



for X < 0, 
for X > 0, 



S{x) = < 



6^/3 (sinh y/—rj — ry)2/3 

1 — cos y/rj 
, Q^^^{^/f] — sinh ^JriYI^ 



for X < 0, 
for X > 0, 



and ^(0) = (3/4) V3. 

The radial and azimuthal Hubble rates are defined by 

rr / X R'{t,r) , , , R{t,r) 
H\\{t,r) := ; and H^{t,r) ■- ^' ' 



R'{t,r) 



R{t,r) 



(A3) 
(A4) 

(A5) 



(A6) 



(AT) 



(A8) 



(A9) 
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By defining the density-parameter function as 



2M(r) 



we can rewrite Eq. flA3p in a form similar to tfie Friedmann equation: 



Hl{t,r) 



HUr) 



( Rojr) 
\R{t,r) 



Rojr) 
R{t,r] 



wliere functions with subscripts correspond to present values, H±o{r) 
Ro{r) = R(to,r). By integrating ( lAllI) . we have 



H±o{r) 



dy 



(AlO) 

(All) 
Hj_(to,r) and 

(A12) 



where y = R/Rq. 

The CR model has the uniform big-bang time tB{r) = 0, and the gauge condition is chosen 
so that R{to, r) = r. In the uniform big-bang model, one functional degree of freedom to 
specify the model remains. In the CR model, this degree of freedom is fixed so that the 
density-parameter function is given by 



(out) 
M 



/Q{0Ut) 

1"m 



(A13) 



where Q 



(out) 
M 



0.7, Q 



(in) 
M 



0.242 and a = 6Gpc. The Hubble constant at the center is 



^0 = H±o{r = 0) = 74kms-iMpc-^ 



Appendix B: transfer function on the CR model 

We use the fitting formula for the matter transfer function T(k) given by Eisenstein and 
Hu 53|. The fitting function defined by equation (16) in Ref. [53| is determined by the four 



parameters as 

T{k) = T{k;n^,n,,h,e2.7), (bi) 

where f2b and Qc are the cosmological density parameter of baryons and cold dark matter, h 
is defined a.s h = iJo/(100kms~^Mpc~^) and the CMB temperature is written as 2.7Q2.7K. 
In this paper, we adopt the following values for these parameters, 

rib = 0.042, fi, = 0.20, h = 0.74, 62.7 = 1.0 , (B2) 

where Qb, and h are chosen to be the same at the center of the CR model and B2.7 = 1.0 
is assumed. This transfer function is plotted as a function of k in Fig. [HI 
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